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Abstract—The effects of phased array antenna patterns on 
copolar correlation coefficient of polarimetric radar are analyzed. 
Analytic expressions are derived and used to quantify biases in the 
copolar correlation coefficient caused by mismatch of beam widths 
for horizontal and vertical polarizations. Misalignment in the 
beams’ pointing direction is also considered. Measurements are 
made on an 8x8 panel of a polarimetric phased array radar 
(PPAR) under development by the Advanced Radar Research 
Center at the University of Oklahoma in collaboration with 
NOAA's National Severe Storms Laboratory. The all-digital radar 
is called Horus and some measurements on a partially populated 
antenna are presented. 


Index Terms—phased array radar, polarimetric radar, all- 
digital radar, weather radar, antenna patterns 


I. INTRODUCTION 


HASED Array Radar (PAR) may be a last step 

improvement in radar technology for observations of 

weather, both to serve scientific investigations and 
operational applications. It is considered for future weather 
radars because it can provide adaptive (i.e., where and when it 
matters) high-quality meteorological observations and thus 
advance understanding of atmospheric phenomena [1] and 
improve warnings and nowcasts. Upgrade of the National 
Network of weather surveillance radars (WSR-88D) to dual 
polarization caused a significant step increase in quality of 
quantitative precipitation estimation as well as capability to 
identify a wide variety of sources causing radar returns [2]. 
Therefore, dual-polarization capability has become a de facto 
requirement for modern weather radars including the PARs. 
Widely accepted polarization mode consists of transmission 
and reception of electromagnetic waves in both Horizontal (H) 
and Vertical (V) polarizations. In the Simultaneous Transmit 
Simultaneous Receive (STSR) mode, adopted by most 
operational and research weather radars, the H- and V-polarized 
waves are transmitted simultaneously thus creating an 
elliptically polarized wave. At the receiving antenna the 
returned waves are separated into channels for horizontal and 
vertical polarizations. Producing accurate polarimetric weather 
measurements with PARs in the STSR mode is challenging due 
to biases induced by the co- and cross-polar antenna patterns 
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when the beam is pointed off the broadside [3] [4]. 

The magnitude of the correlation coefficient between the 
horizontally and vertically polarized returns, hv is a key 
parameter defining the quality of polarimetric radar 
measurements [2], [5]. The standard errors of the estimates of 
polarimetric variables are significantly reduced if pny of the 
weather signals exceeds 0.99 [2] — a basic requirement for 
polarimetric weather radars. Additionally, accurate 
measurements of pnv are crucial for polarimetric detection of the 
melting layer and determination of its height [6] and for 
identification of the areas of hail and quantification of its size 
[7]. Therefore, the requirements for the phy measurements in the 
design of radars for weather observations are strict and 
important. 

The S-Band mobile polarimetric all-digital Horus radar 
system was developed by the Advanced Radar Research Center 
(ARRC) at the University of Oklahoma (OU). Its purpose is to 
rapidly scan storm volumes while yielding high quality 
polarimetric observations. By design the Horus radar system is 
scalable and modular - thus suitable for initial deployments with 
a partially assembled antenna. The companion article 
describing the Horus system architecture is in this journal [8]. 
Considering that the polarimetric performance of a PAR is 
largely dictated by its far-field copolar and cross-polar radiation 
patterns, herein we analyze the patterns of one 8x8 element 
panel of the Horus antenna. Specifically, we investigate the 
impact of differences in the copolar H/V beams on the 
correlation coefficient (pnv) estimates and evaluate these using 
measured patterns. An approximate analysis of the cross-polar 
patterns on hv is in [9]. Furthermore, we extrapolate the 
analysis to a full size (5x5 panels) Horus and a hypothetical 
16x16 panel radar having a beamwidth of about 1° as on the 
WSR-88Ds. 

The analysis of the beam effects on pny assumes a uniform 
distribution of scatterers sizes in regions larger than the 
resolution volume (defined by the beam cross section and range 
depth of the sampling volume). Beams illuminating these 
scatterers cause quantifiable bias in pnv useful for relative 
comparisons. This is unlike the effects on the returned powers 
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(at H and V polarization), differential reflectivity or linear 
depolarization ratio. These, for example, are not biased by 
beams having equal cross section but different pointing 
direction. Bias by beams with different cross sections and 
corresponding gains can be accounted for through calibration 
and measurements (cross sections and gains). This is one reason 
we refrain from considering these biases here, the other is the 
concomitant increase in the article’s length with trivial 
discourse. It is much harder to correct biases in polarimetric 
variables if the beam is pointed out of the array’s cardinal 
planes [10]. 


II. THE ALL-DIGITAL HORUS PAR 


An all-digital phased array architecture is the most advanced 
and flexible PAR technology for enabling enhanced weather 
surveillance strategies envisioned to improve the temporal 
resolution and quality warnings on the basis of weather radar 
observations [9]. Important reasons (e.g., minimal attenuation, 
clear-air detectability) favor an S-band frequency (same as on 
WSR-88D). It can serve a wide community in atmospheric 
science if its antenna is large to form a sufficiently narrow 
beam. The growth of the wireless communications industry 
revolutionized electronic circuits enabling digitization (and 
control) at every element of an S-band array antenna [11]. 
NOAA is exploring digital radar technology [12] and had 
funded the ARRC to develop an all-digital PAR demonstrator 
[8]. 


Fig. 1. The fully digital Horus mobile radar, with 25 antenna 
panels (8x8 dual-pol elements each). More details can be found 
in the Horus System article published in this journal [8]. 


The mobile Horus radar (Fig. 1) is a small S-band 
engineering demonstrator being built by ARRC. Significant 
aim of this project is to mitigate risks (i.e., maturity of the 
technology, maintainability, failure rates of components, 
quality of radar variables, feasibility/tradeoffs for rapid scans 
and adaptive scans) associated with all-digital phased array 
architectures and demonstrate its advanced capabilities. The 
radar has a 4.5° beamwidth suitable for showing functionality 
and overall scalability as well as limited atmospheric research 
about weather close to the radar. Element-level waveform and 
polarization flexibility, imaging, control/calibration methods, 
and power levels are compatible by design with a full-scale 


system for weather surveillance similar to the WSR88D. 

The all-digital array panel defined as an 8x8 mechanical and 
electrical subarray, is the primary building block. The panel's 
tube-based frame provides mechanical support as well as 
conduit for liquid cooling. The “OctoBlades,” eight per panel, 
house all electronics and slide into the panel frame to connect: 
a) the liquid cooling system, b) the backplane, and c) the 
antenna array itself through blind-mate connectors that pass- 
through cutouts in the backplane. The OctoBlades are the heart 
of the overall system implementation [13]. 

A centralized module on each panel is the SuperBlade, 
supporting the active electronic components. This module also 
serves the panel's analog backplane by distributing DC power, 
local oscillator signal, and a reference clock to each of the 
OctoBlades. The radar’s wideband waveforms are switchable 
from dwell-to-dwell in the 2.7-3.1 GHz frequency band. The 
radiators are aperture-coupled, stacked-patch elements that 
provide high-quality polarimetric isolation (better than -45 dB 
cross-polarization suppression after calibration). By digitizing 
at the element level, we have shown that polarimetric 
calibration of the radar hardware can be robustly accomplished 
(see companion Horus article on calibration, [14]). 


II. EFFECT OF H/V BEAM MISMATCH ON (uy 


Herein we present the effects of different beamwidth and 
pointing direction of the H and V beam on the copolar 
correlation coefficient (ny). We ignore contribution by 
sidelobes as well as by cross-polar patterns. To arrive at a 
solution, we make several assumptions and approximations 
(Appendix A) on the basis of which an objective comparison of 
effects is possible. Assumptions are: Uniformly distributed 
scatterers over a volume larger that the radar resolution volume 
Vs [15] (see also the Appendix). Scatterers from Ve contribute 
most to the returned signals. Range extent of Veis same for both 
polarizations. This is expected from matched H and V channels. 
We consider patterns with elliptical cross sections that may 
differ in widths and pointing directions. Therefore, we have 
defined the widths as follows (Appendix). The two-way width 
of the pattern for horizontal polarization along the arc (angle) y 
iS Ony; Om, is the width along the vertical arc (angle) y; for 
vertical polarization the corresponding widths are ov, and oy). 
The two-way width (œ) of a Gaussian pattern is related to the 
one-way 3 dB beamwidth @, via 6, =4,fIn(2)o, [15]. Identical 


relation holds for the vertical (meridional) width and horizontal 
(lateral) width. Because the beams are extremely narrow, the 
overlapping and significant parts of the beam cross section can 
be assumed to be almost flat and parallel to the polarization 
plane (Appendix). In this coordinate system we denote with y 
and y the perpendicular angles in the beam cross section plane 
(Appendix Fig. Al) and integrate the relevant patterns to 
compute the decorrelation by the beam effects. 

The angular (solid angle) cross sections of two overlapping 
beams with different widths and pointing directions are in 
Fig. 2. 
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Fig. 2. Angular cross section of two overlapping elliptical 
beams whose centers are offset by A. 


For Gaussian patterns f and E , the expression for the 
correlation coefficient (Appendix) is 
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where the scattering coefficients at H and V polarizations are 
Shh, Svv, and ( ) denotes ensemble average. In (1) the intrinsic 


correlation coefficient, pnvi is the first term, that is 


(Susu) Q) 
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and the second term is the patterns’ induced bias. The ratio of 
the biased pny (left side of (1)) to the intrinsic correlation 
coefficient (2) defines the fractional bias č caused by antenna 
pattern mismatch 
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In this definition ¢ = 1 signifies unbiased estimate. 

Let the offset of H polarization pattern from zero along the y 
angle be - A /2 and that of V polarization pattern be A/2 (Fig. 
2). Then the equations for one-way patterns are 
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The integrals of (5) as stipulated in (4) are over y and y and for 
the case of interest (narrow beams) can be extended from -œ to 
+o in both dimensions (Appendix). Take circular cross sections 
of different size for the H/V polarizations so THT =a, 


and O, =O, =O,» then integrate (4) to obtain the closed 
form solution for č 
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By inspection we see that with zero offset and same beam 
widths č = 1. Hence there is no bias and the lobes are perfectly 
matched. 

For elliptical beams the f and F are 
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where Ai, A2 are the offsets in the y and y directions. Use of (7) 
and (8) in (1) yields the bias 


For the WSR-88D the pny bias is specified at less than 0.006. 
Because this bias contains the effects of the whole system, we 
choose one half (0.003) to illustrate what a desirable value for 
a PAR antenna should be. Consider circular beam cross sections 
and geometry as in Fig. 2. Then the fractional bias should 
satisfy 0.997 < &< 1.003. By including this inequality in (6) we 
can determine the conditions it imposes on the multiplying and 
exponential terms. Note that the multiplying term has a 
maximum of 1 if on = oy. Further, the strictest condition on the 
fractional bias is if it is closest to 1 (where 1 means no bias). 
This and because the exponential term is 1 if A = 0 means the 
worst (most demanding) condition for the multiplying term is 


0.997 < 2% <1, 
h v 

Clearly the upper condition for the exponential term is also 1 as 

it cannot be any larger. Therefore, the worst (most demanding) 


lower limit for this term is 


(10) 
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0.997 < exp - (11) 


A? 

A(o; +o) 
One can see that in this case the condition (10) depends on owov 
and (11) depends on one more parameter, A/o. For elliptic 
beams the following four parameters Ohy/ Owy, Ai/Ovy, Ony/ Ov, 
and A2/Ohy determine the bias in Pay. 

Next for simplicity and illustration assume matched pointing 
directions so that A = 0 in (6) and the mismatch in beam widths 


1S 
o, =(lté)o, (12) 


where € is a number much smaller than 1. Inserting (12) into (6) 
one obtains the following relation, 


+e)o? E 
aaoi EE e og 
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Assuming beam mismatch is the only cause of pav bias we need 
to find the largest e for which the inequality (13) holds. The 
solution for (1e) is 


2 
dsos AE S a4) 
consequently 
2 
pe N = =i, (15) 


As example assume a drop in correlation from 1 to 0.997, 
hence € = 0.997 and (rounded) €12 = +0.08. Then if the 3 dB 
beamwidth (H polarization) On = 1°, the @, would be 
constrained to between 0.92 and 1.08°. This can be also cast as 
difference at an L (dB) level below the peak of the pattern as 
follows. Note that the angle with respect to beam center at L 


(dB) below the peak is 
O N2L o, V2L 


4,Jin2[10log(e)] 10 log(e) | 


At mismatched beams (10) holds and we can compute the 
difference 520 (dB) between matched (e = 0) and not matched 
(e # 0) patterns at the L level as 


1 L 1 
(+e (1-67 


Thus if 1+ = 1.08 and L = 20 dB, the 629 = 2.85 dB. If hv = 
0.999, 520 = 1.7 dB so it is fair to state that at the 20 dB level 
the difference in the patterns should be less than 2 dB, provided 
that at the peak the gains have been matched. 

In a practical array it may not be possible to have a perfectly 
circular beam cross section but an approximately elliptic one, 
especially in the case of a phased array radar pointing off 
broadside. Let us look at elliptic Gaussian patterns with widths 
Ohy, Ohy and Ow, and Ow as explained earlier and depicted in 
(Fig.2), but with both H and V lobes centered at the origin. Then 
the fractional bias č becomes (Appendix) 
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Next consider a special case where the lengths of the longer 
axis are equal and also the lengths of the shorter axis are equal. 
That is Oh = Ov and Ow = Ohy. These are implied by the 
symmetry of the patch’s shape. Moreover, if centers of the two 
cross sections coincide then rotation by 90° of one of these 
would produces the other one and (18) reduces to 
ee 2 Oy Fy (19) 

(o, + o) 

This equation represents the estimated bias by antenna patterns 
with elliptical beam cross-sections. By setting A = 0 in (6) we 
see that (19) is same as bias by circular beams with different 
widths. The reason that the two geometries yield the same result 
is in the underlying assumptions about the scatterers. These are 
uniformly distributed over a volume much larger than the radar 
resolution volume. 

In the next section, we use measured patterns of an 8x8 active 
Horus antenna panel to estimate their influence on phy 
estimates. 


IV. EVALUATION OF HORUS ANTENNA PANEL PATTERNS 


One antenna panel of the scalable Horus array was fully 
populated with electronics that make up the transmit and 
receive signal paths and was assembled in a near-field chamber 
for testing (Fig. 3). This setup enables the panel to operate as a 
radar (albeit with a small aperture) so that we can test its 
performance in the STSR mode for dual-polarization 
measurements. The near-field scanner is comprised of two 
motorized Velmex BiSlide assemblies, one Velmex VXM 
Stepper Motor Controller, an S-band open-ended rectangular 
waveguide probe (OEWP), a Newport “optical breadboard” 
base, and RF absorber. These features allow measurement of 
antenna patterns in the H and V polarizations. Two-axis motion 
is achieved by controlling one stepper motor per axis. The 
motor controller takes in formatted commands and drives the 
stepper motors to achieve the desired movements. The S-band 
OEWP is mounted to the near-field scanner. The probe is 
covered with a collar of RF absorber. 


Fig. 3. Near-field scanner setup to fully characterize the H and 
V antenna patterns of the Horus panel. 
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Fig. 4. Normalized (left) H and (right) V far-field Horus panel patterns derived using the near-field to far-field transformation. 


The near-field scanner is designed so that the system’s (0,0) 
position is in the exact geometric center of the array face. 

Near-field calibration is accomplished using the so-called 
park-and-probe technique to measure amplitude and phase of 
fields radiated by each antenna element, generate and apply 
alignment weights, and then verify the resulting alignment. A 
full set of transmit patterns requires four separate data 
collections, one for each combination of array polarization and 
OEWP orientation (0° or 90°). Although the panel hardware is 
capable of dual-polarized transmit operation, the OEWP is 
linearly polarized. Therefore, it is impossible to separate the H- 
from the V transmit response if transmitting both polarizations. 
Receive patterns can be collected with only two data sets. The 
panel hardware can receive both polarizations simultaneously 
and feedback the data separately. Multiple beam angles can be 
collected simultaneously in the receiver. Considering that a 
thorough alignment procedure was conducted and resulting 
H/V beam alignment was verified, hereafter, we neglect 
misalignment and consider only the beam shape mismatch. 
Beam shapes are determined by the design of the microstrip 
patches exciting each polarization and are more challenging to 
correct after the array is designed and manufactured (in contrast 
to alignment). A deeper investigation on the impact of H/V 
beam misalignment on polarimetric PAR measurements is in 
[16]. 

Determining the far-field pattern of an antenna from near- 
field measurements requires a mathematical transformation and 
correction for the characteristics of the measuring antenna (i.e., 
the probe). We apply the near-field to far-field transformation 
via mathematical expressions derived by [17] to estimate the 
antenna far fields, similar to the procedure in [14]. 

Field transformations expressions are applied to the Horus 
near-field measurements to predict the far-field antenna 
patterns of this 8x8 array (i.e., one Horus panel), as in [18]. 
Normalized co-polar H and V far-field patterns are presented in 
Fig.4. Dotted contours on the copolar H and V patterns indicate 
the half-power beamwidth (-3 dB). 


A comparison of these patterns indicates excellent mainlobe 
agreement between the H and V polarizations up to about -20 
dB. The sidelobe structure for each polarization appears to be 
symmetric about the mainlobe for the horizontal and vertical 
cuts. A small difference (<0.7 dB) is observed between peak 
sidelobe levels in the H and V patterns. This may be caused by 
finite array effects such as diffraction and mutual coupling [19]. 
Similarly, the diagonal plane cuts show well matched H and V 
mainlobes (up to about -20 dB), although sidelobe levels are 
lower (as expected) due to the rectangular antenna geometry. 

To evaluate the impact of these far-field patterns on pny 
(eqs.12 to 17), the beamwidths from cross sections at -3, -10, 
and -20 dB levels are extracted from the measurements and 
presented in Fig. 5. The widths measured in y and y at the -3, - 
10, and -20-dB levels are summarized in Table I. 

Next, we compute the fractional bias č that the panel’s 
mainlobe would inflict to the pny estimates. We can use (18) for 
that purpose by substituting the values from Table 1 because the 
Gaussian model allows use of widths from an arbitrary level 
below the mainlobe peak. Our results from the three levels 
(Table 1, last row) for all practical purposes validate this 
assumption. The slight difference between č computed from the 
width at the — 3dB level and the widths at the -10 dB or -20 dB 
level are negligible. Moreover, the difference €(-3) - &(-L) 
increases with the increase in the level’s magnitude L. This is 
likely due to the increase in deviation of the mainlobes shapes 
from the Gaussian at larger levels below the peak. Therefore, it 
is appropriate to take the widths at the -20 dB level as these 
impose more stringent requirement. This also follows from the 
fact that statistical errors in pnv are not affected if SNR is larger 
than 20 dB [2]. 

Now we estimate the pny bias a PAR may cause if it has a 1° 
beamwidth as on the WSR-88D [15]. For that we take the 
widths at the -20 dB level and assume that the same absolute 
deviation from the mean width as in the one panel case would 
occur. The mean for the panel from table 1 forth column is 
25.681° and the mean width for the 1° beam width antenna at 
the -20 dB level from (16) is 1.825°. Then the ¢(-20) of the 1° 
antenna from (18) is 0.9874 corresponding to a bias of 0.0126. 
For this is out of specification value to occur all panels on the 
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antenna would have to cause identical bias that would add 
“coherently”. Itis more likely that the panels would differ and 
the deviation of the widths from the mean would be random. 
Let us consider that case and assume the antenna consists of 
16x16 panels. Then the variations of the widths would be 
approximately 4 times smaller than the deviations 
corresponding to the widths at the -20 dB level in Table 1. The 
resulting č = 0.9993 and the corresponding bias of 0.0007 is 
negligible. We do not know if coherent or random addition of 
widths fluctuations will apply to the fully populated Horus 
antenna. Nonetheless, we plan to determine the bias Horus 
mainlobes impose once the antenna is assembled and pattern 
measurements are made in the ARRC’s recently installed near- 
field scanner [14]. It may be possible to extrapolate these 
measurements to ones expected on a 1° beamwidth Horus type 
antenna. 


TABLE I — Measured Pattern Widths (Boresight) 
Pol. \ Level -3 dB -10 dB -20 dB 


H %) IDA Sia 20.964° 25.662° 
H (y) 12.558° 20.802° 25.42° 
V 12.556° 21.052° 25.716° 
V) 12.581° 21.112° 25.928° 
Bias ¢ 1.0 0.99998 0.99996 


The results derived here are obtained using a single Horus 
panel. From the expressions derived in the previous section and 
the measured widths, we estimate that the antenna will 
introduce at least an order of magnitude lower bias in hv than 
the requirement for polarimetric observations. Although a 
larger aperture would result in a narrower beam, it is expected 
that the H/V beam shapes will be preserved providing that 
accurate alignment calibration is performed. 
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Fig. 5. Beam cross-sections at the -20, -10, and -3-dB level 
below the peak. 
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Fig. 6. Beam cross-sections at -3 dB as a function of 
steering angle, derived from NF measurements obtained 
using the setup presented in Fig. 3. The beams are 
steered in elevation, going from -42° to +42°. Dashed 
lines correspond to H and dotted lines to V. 
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V. PATTERN PERFORMANCE OFF BORESIGHT 


It is known that co- and cross-polar antenna patterns of 
polarimetric PAR induce biases in the polarimetric variables, 
which increase as the beam is steered away from the broadside 
[20] [21]. Here in subsection A we present panel patterns off 
broadside and associated bias in pny. In B we consider a full-size 
Horus antenna and the bias it induces in pny as function of the 
pointing direction, and in C we present data obtained with an 
1x5 panels Horus antenna. 


A. Panel patterns off broadside 

We consider Horus panel patterns electronically steered along 
the vertical principal plane. Specifically, a set of 15 independent 
measurements were taken per polarization as the beam is 
steered from -42° to 42° at 6° steps in elevation. In these 
measurements, the S-band probe transmits, and the Horus array 
receives. First, the probe was aligned to transmit horizontal 
polarization, then, to transmit vertical polarization. Co- and 
cross-polar antenna patterns are collected for each polarization 
and steering angle. 

For each measurement, the beam cross-section at the half- 
power beamwidth is extracted and the pattern widths in y and y 
are measured. Measured beam cross-sections are shown in Fig. 
6. The H and V cross-sections are represented with dashed and 
dotted lines. As expected, when the beam is steered away from 
the broadside, the pattern width increases approximately as 
1/cos(@), were @ is the steering angle. At broadside, the widths 
are ~12.6°, and at the elevation steering angles of +42°, widths 
are ~16°. 

Qualitative comparison of the H and V cross-sections 
indicates that widths are well matched up to about +30°. 
Beyond these steering angles there seems to be a small pointing 
offset between H and V, the worst one being < 0.1°. 
Nevertheless, beam shapes appear to change similarly in 
general, which means that biases induced by the copolar 
patterns in hnv are small. Beam widths measured at the -3 dB 
level are shown in Fig. 7. 

It is apparent that the V pattern beamwidth increases more 
than the other patterns along y, by about 1° at the +42° steering 
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Fig. 8. Impact from copolar Horus panel antenna 


patterns on Phy derived using (18) and the measured 
widths shown in Fig. 7. 


angles (Figs. 6 and 7). We hypothesize that the reason is that 
the V patterns are aligned with the steering plane (vertical), 
which in a small array (8x8) may cause larger differences in the 
pattern shapes due to diffraction effects [22]. These measured 
beamwidths are used in (18) to estimate copolar antenna biases 
introduced by the Horus antenna in pny, shown in Fig. 8. Results 
indicate that biases are smaller than 0.0005) within the 
elevation (0) span of + 20°. They increase slightly as the beam 
is steered to +42° where the maximum of ~0.0015 occurs. This 
is about half of the desired pny bias stated earlier (0.003). 
Measurements in the azimuth plane yielded very similar results 
hence are not presented. 


B. Horus antenna patterns off broadside 

Although panel patterns were only measured in the elevation 
plane, the intent is to evaluate the performance of Horus PAR 
in a larger scan sector, steering the beam in two dimensions. 
Further, it is of interest to evaluate the performance of the full 
Horus array (4x4 active panels, 32x32 elements). Herein, we 
use element pattern measurements obtained in the ARRC’s 
precision anechoic chamber and in the spherical far-field scan 
mode. From these we quantify biases of the full Horus array as 
the beam is steered in azimuth +45° and in elevation from 0° to 
20°. A simulation following the approach proposed by [10] [16] 
is used. It combines the effects of measured radiation patterns 
with simulated time-series signals, applying the well- 
established backscattering matrix model [23]. The model 
includes bulk statistical properties of scatterers within a 
resolution volume coupled with electromagnetic wave 
propagation and radar system effects. Effects from both the 
copolar and cross-polar antenna radiation patterns are 
considered. Simulation results are shown in Fig. 9. 

It can be seen that biases induced by the Horus antenna in pny 
agree with those in Fig. 8 (along the vertical principal plane, 
where cross-polar patterns are < -50 dB), and are below the 
requirement [24] (i.e., bias in pny < 0.006) through the entire 
scan sector. This is due to the well-matched copolar beam 
patterns and the ultra-low cross-polarization levels of the Horus 
antenna [25] (< -50 dB in the principal planes of both the full 
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Fig. 7. Beamwidths measured from the cross-sections in 
Fig. 6. Note that the curves for ony, and Ow are partially 
overlapping. 
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antenna and the panel). The integrated main lobe power to cross 
polar pattern power measured on the panel is 47.9 dB. 

In summary, our quantitative analysis based on measured 
antenna patterns suggests that any adverse impact of the Horus 
antenna on ppv estimates will be negligible. We expect that the 
shapes of H and V beams are preserved for the full Horus array 
to produce high-quality polarimetric weather measurements. 


Bias in Pry 


20 0.005 
0.004 
15 
0.003 
10 
0.002 
5 
0.001 
0 0 
-40 -20 0 20 40 


Azimuth (°) 


Fig. 9. Antenna-induced biases in pny as a function of 
steering angle, derived from measured element patterns for 
the active part of the full Horus array (4x4 panels). Note that 
biases induced are below the requirement (i.e., bias in pny < 
0.006) through the scan sector of interest. 
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C. Evaluation of Horus pw Estimates 

On the evening of 16 February 2023, the Horus radar was 
deployed near the Radar Innovations Laboratory in Norman, 
OK. Soon after, at approximately 1:14:22 Z, a quasi-linear 
convective system of storms was moving South-West to North- 
East passing south of Norman. Severe thunderstorm, flash 
flood, and tornado warnings were issued by National Weather 
Service (NWS) forecasters throughout the event!?. We 
collected polarimetric data from this event with the Horus radar 
for approximately 30 minutes. 

Five Horus antenna panels (5x1) with corresponding 
electronic components (i.e., transmit/receive chains) were 
assembled into a vertical column of the Horus array structure at 
that time. We performed polarimetric calibration using the 
mutual coupling technique [26] and collected meteorological 
data for engineering evaluation. We scanned the beam in the 
vertical principal plane through a convective storm cell, going 
from 0° to 30° in 0.5° increments and obtained field of the radar 
variables. A pulse-repetition time (PRT) of 1 ms with 128 
samples per dwell was used, resulting in a scan time of ~8 s. 
Pulse compression waveforms with non-linear frequency 
modulation were used, with a pulse width of 80 us and a 
bandwidth of 5 MHz. Note that the progressive pulse 
compression technique was implemented to mitigate the pulse- 
compression blind range [27] [28]. Range-time samples were 
produced at a rate of 15.625 MSPS, resulting in a range 
sampling interval of approximately 10 m. 


' See: https://www.spc.noaa.gov/exper/archive/event.php ?date=20230215 


Pertinent to this paper are the signal-to-noise ratio (SNR), the 
correlation coefficient and the differential phase (Ppp) because 
it is the argument of the correlation coefficient, which are 
shown in Fig. 10. A qualitative analysis of the fields shows SNR 
values going up to approximately 40 dB at close range. 
Generally smooth fields of estimates were produced, as 
expected from volumetric scatterers encountered in 
precipitation systems that generally exhibit spatial continuity. 
Changes in ®pp are consistent from those in dish-based weather 
radars, whereby larger increases are observed in regions with 
high reflectivity (i.e., proportional to SNR at constant ranges). 
This is caused by the change in phase encountered by H and V 
waves as they traverse a precipitation core, where raindrops 
tend to be larger in the horizontal dimension (due to gravity and 
drag force), and therefore the H wave propagates slower than V 
(i.e., larger phase change) producing Ppp increases along range. 
The field of pny appears to have most values close to 1, 
indicating most scatterers are of the same type (i.e., liquid 
raindrops or frozen hydrometeors). Regions with lower pny - 
values may indicate mixed precipitation (e.g., water and ice), or 
a more general mixture of particles (e.g., debris). 

To quantify polarimetric Horus performance, we produced 
histograms of pnv values and of differences between fields of ppv 
derived from consecutive scans, presented in Fig. 11. To 
mitigate estimation biases, only data with SNR = 20 dB are 
selected. Histograms are normalized to approximate a 
probability density function. The height of each bar represents 
the ratio of the number of observations in bin divided by the 
product of the total number of observations by the width of a 
bin (in this case 0.0005). The area of each bar is the relative 
number of observations, and the sum of the bar areas is 
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Fig. 10. Polarimetric weather measurements from a 5x1 

panel Horus array (40x8 elements). Data were collected at 

approximately 1:14:22 Z pointing at 190° azimuth. Panels 

are (top) SNR, (middle) ®pp, and (bottom) pnv. 


? https://nwschat.weather. gov/lsr/#/202302152220/20230216083 1/0100 
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less than or equal to 1. The top histogram shows that most pnv 
vales are concentrated between 0.99 and 1. This indicates H and 
V beams are well matched in shape, therefore antenna-induced 
biases are negligible. This is consistent with our analysis 
presented in Section IV. Furthermore, a zero-mean histogram is 
produced using differences between fields of pny from 
consecutive scans. The distribution appears to be a bilateral 
exponential or (Laplace distribution), with a standard deviation 
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Fig. 11. Histogram of (top) pnv values, and (bottom) 
differences between two consecutive fields of pnv. Note that 
an SNR treshold of 20 dB is used to select data used in these 
histograms. 


of 0.0085. This indicates stability of the system. Although 
standard errors in pnv are mostly dependent of SNR, dwell time, 
and spectrum width, these histograms show that the Horus radar 
can produce polarimetric estimates of high quality. 


VI. CONCLUSION 


This article presents progress towards evaluating the 
performance of the Horus radar for polarimetric weather 
observations. Horus is a digital, polarimetric, S-band phased 
array radar being developed at the ARRC with support from 
NOAA. The Horus antenna aperture consists of a 5x5 
arrangement of panels (of which not all are active), each with 


8x8 antenna elements. With approximately 1024 active 
elements in total (4x4 panels), this proof-of-concept all digital 
system will have a 3-dB beamwidth of 4.5°. The Horus radar 
system has been integrated and initial deployments for 
polarimetric weather measurements are ongoing. 

A theoretical evaluation of the mismatch between H and V 
beams was presented, with a focus on its impact on the copolar 

correlation coefficient (Pnv). Two mismatch parameters were 
considered: a pointing offset between beams and the ratio 
between the 3-dB widths in x and y (for H and V). Assuming 
magnitude and phase alignment is performed, and relatively 
narrowband waveforms are used (i.e., no beam squinting), the 
pointing offset would be mitigated. However, the widths of the 
H and V beams are determined by the element patterns and 
may be more challenging to adjust. Therefore, we derived 
theoretical expressions assuming no pointing offset but 
potentially different H/V beam widths. 

In our measurements deviations between the widths at -20 
dB below the peak are larger than at higher levels and are 
easier to estimate that deviations at the — 3dB level. Therefore, 
and because pn is not affected at SNRs larger than 20 dB, the 
widths at -20 dB from the mainlobe peak approximate better 
the effects on bias. Using these widths, we estimate that a 
complete Horus antenna would cause insignificant bias. We 
extrapolated the results to a Horus type PAR with a beamwidth 
of 1° like on the WSR-88D. The result indicates that if the 
deviations of single panel pattern widths sum coherently the 
bias would be unacceptable. Nevertheless, if these biases are 
incoherent then the resulting bias by the 1° antenna would 
satisfy the requirements. Incoherent bias means that the 
variations of the widths from panel to panel are independent 
which is more likely. Otherwise, all panels would have 
identical “flaws” and alignment error of panels on the antenna 
structure would be same. 

Using near-field measurements from a fully populated 
active Horus antenna panel, we derived cuts of the antenna 
patterns in the x and y directions and for both H and V 
polarizations. A qualitative comparison of the H and V patterns 
indicates there is a high level of similarity. Moreover, the 
measure integrated copolar power of the panel to the cross 
polar power at broadside and within the main lobe is 47.9 dB. 
This and our quantitative analysis of the two (H, V) main lobe 
patterns, based on their measured widths, suggests that the 
impact of the Horus antenna on hnv estimates will be 
negligible. We expect that the shapes of H and V beams are 

preserved for the full Horus array and thus will produce high- 
quality polarimetric weather measurements. 

We collected polarimetric weather data with the 5x1-panel 
Horus array and evaluated its performance on the estimation of 
Pov. Relatively smooth fields of polarimetric weather estimates 
were produced (as expected from a precipitation system), from 
RHI scans collected over a period of ~30 min. The majority of 
Pnv estimates are between 0.98 and 1, i.e., 91% of estimates with 
SNR = 20 dB. This is expected from hydrometeors of the same 
type and it reflects on the high level of copolar H/V beam 
matching attained by the antenna. Furthermore, since these data 
were collected at different electronic steering angles along the 
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vertical plane, this shows that beams are also well matched 
along the vertical cardinal plane. These observations are 
consistent with the results obtained in Section IV. Furthermore, 
element-level mutual coupling measurements are being 
recorded regularly and will be archived. In the future, we plan 
to analyze these measurements to quantify magnitude/phase 
drifts as a function of time, which impact polarimetric 
calibration and are a function of hardware stability. 
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APPENDIX 


Model for antenna patterns 


We present a succinct derivation on the basis of which the 
effects of mainlobes on the pny are evaluated. For characterizing 
antennas, the Ludwig 3 coordinate system [29] is convenient as 
it is tied to the antenna. Herein we use it as follows. The 
coordinate z we define along the beam center, the vertical 
coordinate we label as v and the horizontal coordinate 
perpendicular to both v and z we label with h (Fig. Al). To 
compute the polarimetric variable, one needs to evaluate 
integrals over the radar resolution volume Ve centered at range 


11 


ro. Ve is defined by the beam cross section and the range 
weighting function [15]. We will use the arc lengths rody and 
rody for deriving various second order moments of voltages 
representing backscattered electric fields. 
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Fig. A1. Coordinate system used in defining various patterns 
and orientation of scatterers for deriving bias in pnv. The view 
is from above the polarization plane (h, v). The beam center is 
along the coordinate z which is also the direction of 
propagation. The horizontal plane is defined with z, h. The 
vertical axis is v and the vertical plane is defined with v, h. The 
projection of a radial (not drawn) defining a point within the 
beam on the z, h plane is r,, and r, is a similar projection but on 
the z, v plane. 


The angle y is between the projection r, of the radial defining a 
point within the beam onto the v, z plane and the z axis [Fig. 
Al]; the angle y is between the projection of r, onto the h, z 
plane and the z axis. The elemental arcs we defined are 
orthogonal at the beam’s center and almost orthogonal in the 
vicinity of the center. Consequently, we can assume 
orthogonality within pencil beams. The standard orthogonal 
system in Ludwig 3 consists of angle ¢ on a circle in the v, h 
plane centered on the z axis; 0 is the angle between the z axis 
and a radial on the cone (the cone axis is along beam center). In 
that system a standard notation for a Gaussian main lobe one- 
way pattern [15] is 
2 


rop-a -2 ) for 80 €[0,z] and ¢ e[0,27]: (A1) 
O9 


where oy is the width of the two-way pattern, f*(0,ġ) and is 


related to the one way 3-dB width via 8, = 4,/In(2)o, , [15]. 


The integral of (Al) over the Gaussian lobe equals 
10; | (8In 2)and is part of the weather radar equation (see 


[15], eq. 4.13). 
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Next, we demonstrate that the integral of (A1) over the angles 
0, ¢@ can be replaced with the integral over angles y, y of the 


function Fe (X, y) which explicitly is 


f f n= f fel- x E E 


7 
wee oe 201, 


where the first index on f Lae the intended polarization, 
i.e., port for horizontal polarization is exited and horizontal 
polarization (second index) is achieved. This port also 
generates a vertically polarized wave through the cross polar 
pattern. That one we do not consider as our primary goal is to 
quantify the main lobe (copolar) effects. The first index on ony 
indicates the pattern corresponds to the horizontal polarization 
and the second indicates that on, is the pattern width along the y 
angle (arc); similarly, on, is the width of the pattern for 
horizontal polarization but along the almost perpendicular angle 
(arc) y. The integration limits in (A2) have been extended to - 
œ, oo similarly to the limit of the (A1) integral. Assume circular 
pencil beam so ony = ony = oo and integrate (A2) to produce 


= m0; /(8In 2) exactly same is the integral of (A1). We 


will use the expression for f$ (x,y) in (A2) and an analogous 


(A2) 


xdy. 
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one for the vertically polarized wave to compute the correlation 
coefficient. 
We start with the transmitter voltages V; , V’ for the H, 


V polarizations that generate incident electric fields on a 
scatterer located at ri (range ri, and angles yi, yi). To isolate the 
effects solely attributable to the main lobes we consider only 
the copolar scattering coefficients sp, (in the h direction 
coincident with the arc y) and svv (v direction coincident with 
the arc y). We also ignore the cross polar pattern functions Fhv 
and Fw» and express incremental received voltages 
corresponding to a single scatterer as 
t 
e 


ôV, Fin Ho 7 WO) Sm (È) 0 

ov, 0 Fy Lr VW) SO) 

where index i identifies the scatterer and C contains the range 
dependence and other factors in the radar equation so that the 
right and left sides of (A3) are consistent and dimensionally 


correct. The pattern weighting functions for the two 
polarizations are 

Finn EN En try (A4a) 
and 

FY) ENE, GY) bj (A4b) 


Antenna directivities are gn, gv, and fh, fy are normalized to their 
peaks. For compactness we use the following equalities and 
abbreviations 


ROY =ORO=i ae i 
(AS5a) 
and 
2 2 2 2 p2 -x° 40%, -y doy, 
F; NE fa Of, = ffo = e 


(A5b) 
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The received voltages from scatterers principally inside the 
resolution volume are superposition of incremental voltages, 
therefore 


Vi = Cr Fay | wr) | ae (i). 


We compute the correlation coefficient from the ensemble 
averages of the second order moments 


- A . 


By definition the relation between the moment (VAV) and the 


(A6) 


(A7) 


antenna radiation patterns is 


(vi K)e(ER ao EA AeA) LAC) |) Cae SOs) 


(A8) 
The range to the resolution volume center is ro and w(r) is the 
range weighting function including sinusoidal dependence 
within the pulse, i.e., w(r)=|w()|exp(-j4ar/1). The summation 
of |w(r)|? is over the interval ôr centered on ro where |w(r)|? is 
substantial. For comparative purposes we assume scatterers are 
uniformly distributed, have locally invariant statistical 
properties and no differential phase upon scattering. Then 
(sia (s,,@) = (Sasa) and can be pulled out of the sum. After this 


simplification, commute the expectation and summation and 
separate the independent parts to express (A8) as 
(VV) = C (Sin Se) D (Fa ie 4) Fe Ko OIW [LG oP"), 


ij 


(A9) 
Moreover, uniform distribution implies that differences rj- ri 
are also uniformly distributed so that the second ensemble 
average in (A9) differs from zero only if ri= rj (1.e., i=j) and this 
condition reduces the double sum to the single sum 


(Viv) =C (Sin Sw) 2 Fin (io 7) Ea (Z7) | wr) [ 


(A10) 
Replace the sum with integrals (for details see [15], sec 4.4] 
over the resolution volume. The integral over the angles is 
separable from the integral over range; the latter one is a 


constant dependent on wr) , which in turn is a function of the 


pulse shape and receiver impulse response [15] hence has no 
bearing on the correlation. Use the expressions in (A3) but for 
beam cross sections specific to each polarization, and write 


© œ% 


(VEV YE (SmS) n8, | | MAGNA UN xdy = 


œ œ 2 pi 2 y? 
x y x = 
(s Sin’ Sy Yag, Í fe 4o, 4c? J 40? 40 2 A jazar 


—% —0 
(s Shn vv ) 


AIO), Tn OvO vy 


Ehv 
Jae +0, \o, +o, 
(All) 


Via the same reasoning we get the moments in the denominator 
of (A4) yielding 
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5 ane ‘ 
(Iv, F< (sane fem 2s Jef- b- 


—00 


ea ce ) 882A, Ony 


(A12a) 
and 
7 e 2 
(iv, Pho (si) f exp e x | exp T y= 
= (8%) 8n8, 200,0 
(A12b) 


Substituting (A11), (A12.a), and (A12.b) in (A8) the correlation 
coefficient reduces to 


Uns) v2 °° 
(af) (sf) 


| ROR dz S| EOY 


Í fiaz] find i far | fdr 


(i) 24 PnP ny PvP vy 
1/2 1/2 i 
2) 2 2 2 2 2 
(lsm ) ( ) Je, Tiley (or, +O, 


Introducing the compact notation (6a) and (6b) further 


simplifies (A13a) to 
3\ 1/2 32 
(lsm) (lwt) 


Sind 
| ROROA S KOR Mar 


Pry = 
S 


vv 


S 


vv 


(A13a) 


Paw = 
S 


vv 


1/2 * 


Í fe (nar | fiar Í ETA fi dy 


(A13b) 
As explained in the main text the intrinsic prvi is the fraction 
containing the ensemble averages of scattering coefficients 
products. The second fraction with integrals over the pattern 
lobes is the fractional bias č caused by the patterns 


[ROR dz) ROR MAY 
é = 90 —0 


1/2 yar 
jÍ fiaz f siar] Í fiwa | oar] 
(A14) 
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